Practionalization and Fermi surface volume in heavy fermion compounds: the case of 

YbRhs Sis. 
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We establish an effective theory for heavy fermion compounds close to a zero temperature Anti- 
Ferromagnetic (AF) transition. Coming from the heavy Fermi liquid phase across to the AF phase, 
the heavy electron fractionalizes into a light electron, a bosonic spinon and a new excitation: a 
spinless fermionic field. Assuming this field acquires dynamics and dispersion when one integrates 
out the high energy degrees of freedom, we give a scenario for the volume of its Fermi surface 
through the phase diagram. We apply our theory to the special case of YbRh2(Sii_i Gea:)2 where 
we recover, within experimental resolution, several low temperature exponents for transport and 
thermodynamics. 
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In some heavy fermion compounds, pronounced devia- 
tions from the conventional Landau Fermi liquid behav- 
ior have been reported when the compounds are tuned 
through an antiferromagetic (AF) quantum critical point 
(QCP) [1]. Generically, the specific heat coefficient is 
seen to diverge at the QCP, showing at least logarithmic 
increase [2-5] as the temperature is decreased. The re- 
sistivity is quasi-linear in temperature [6-9]. NMR and 
/i-SR studies, as well as neutron scattering measurements 
for one compound [4] show that the spin susceptibility ac- 
quires some anomalous exponent [6,7]. Here we focus on 
the special case of YbRh2(Sii_a: Gex)2 doped with Ge 
to reach a QCP at x = 0.05. In this compound, linear 
resistivity is observed, at the QCP, over three decades of 
temperature, from 10 K to 10 mK [10]. At the same time, 
the specific heat coefficient shows an upturn below 300 
mK from its original logarithmic increase, C/T ~ T~" 
with a ~ 0.33 [7,11]. The entropy associated with the 
upturn is of the order of 5 % of the total spin entropy 
of the Yb atom. Curiously, the upturn in specific heat is 
un-correlated with transport, since one doesn't see any 
kink at 300 mK, or any deviation from the linear slope 
in the resistivity. An NMR study on the pure compound 
[12] shows that the relaxation time on Si 1/TiT - T^^/^ 
when the system is driven close to the QCP with applied 
magnetic field, suggesting that J2qX"{QT^)/^ ~ T~°-^. 
The Griineysen parameter, ratio of the thermal expan- 
sion coefficient and the specific heat is shown to diverge 
with an unusual exponent, i.e. F — (3/C ~ T""'^, com- 
pared to the case of CeNi2Ge2, well understood within a 
Spin Density Wave(SDW) scenario [13]. Latest measure- 
ments of the Hall constant [14] suggest that the Fermi 
surface volume is increasing, when one goes from the 
AF to the field induced heavy Fermi liquid phase. A 
recent study of the heavy Fermi liquid phase shows a 
generic scaling in B /T in the transport and specific heat 
[11]. One also observes that the Kadowaki- Woods ratio 



- K — where A is the coefficient of the resis- 

tivity and 7 — C/T\t^o is the specific heat coefficient- 
increases in approaching the QCP. Likewise the Wilson 
ration - W = xllW^f) where x is the bulk susceptibility- 
shows a dramatic increase. In contrast, the ratio of A/x^ 
stays constant over the whole magnetic field range. We 
deduce from the fact that the A coefficient doesn't feel 
the upturn in specific heat, that the bulk susceptibility x 
doesn't feel it either -their ratio being constant [15]. In 
other words, the excitations responsible for the upturn 
seen in C/T don't couple to the bulk susceptibility. 

In this paper we introduce an effective theory for heavy 
fermion compounds, which, in the particular case of 
YbRh2Si2, accounts for the whole set of experimental 
evidence quoted above. 

Our starting point is the Kondo-Heisenberg lattice 
model, where we use a Schwinger representation for 
the spin of the impurities, which we refer to as the 
bosonic Kondo-Heisenberg (BKH) model. The BKH lat- 
tice Hamiltonian 

H = Hc + Hk+Hh , 
where He ^kfLfka , 

ka 



Hh = Jh E 



(1) 



describes the conduction band (i?c), the Kondo coupling 
between local moments and the conduction electrons at 
site i (Hk), and the super-exchange between neighboring 
spins (Hh)- When we formulate the BKH model as a 
functional integral, we can decouple the fields as follows 
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FIG. 1. Evolution of the Fermi surface volume from the AF 
phase to the heavy Fermi liquid phase. D is the bandwidth. 



FIG. 2. Two possible phase diagrams in heavy fermion 
compounds. 
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where the on-site bond variable in the first term is a 

Grassman field which doesn't carry a spin and the bond 
variable in the second term has been chosen following an 
SP(2) decomposition of the interaction [16]. 

We believe the fermion field % is a good decoupling 
of the Kondo interaction in the sense that it describes 
the set of elementary excitations close to the QCP [17]. 
Through coupling to the spinous ( hka) and the itinerant 
electrons {fka) the x-fermion Kondo bond states acquire 
some dynamics, damping and dispersion. Within this set 
of variables, the formation of the Kondo heavy quasipar- 
ticle is described with the formation of a "band" for the 
X-fermion bound state. To see how this happens it is 
instructive to step back to the impurity case. 

Magnetism in the bosonic Kondo impurity model has 
been investigated in the past [18,19], but here we focus on 
the fate of the fermion chemical potential [20]. From 
(2), one notices that, for the antiferromagnetic Kondo 
effect {Jk > 0), the energy level is negative, thus full, 
and the chemical potential fi = —1/Jk flows to zero at 
strong coupling. The x-fermion becomes massless. This 
singularity accounts for the increase of the Fermi surface 
volume by one. Oppositely, in the ferromagnetic case 
{Jk < 0), the energy level is positive, thus empty, and 
the chemical potential ^ ~ —1/Jk flows to infinity at low 
energies, leaving the Fermi surface volume unchanged. 

In the Kondo lattice, the x-fermion acquires some dy- 
namics, damping and dispersion. The competition be- 
tween AF fluctuations and the Kondo screening affects 
the formation of the x-band. Moreover, one can imag- 
ine that the lattice accommodates partial screening of 
the impurity spin through the x-dispersion: the "un- 
screened" part of the spin -by analogy with ferromag- 
netic Kondo- has positive energies, while the screened 
part of the spin has negative ones. Deep inside the AF 
phase (diagram a in Fig. 1 ), the spin of the impurities are 
partially screened and the x-ferrnion bandwidth is large 
- partially filled band, with a rather "small" Fermi sur- 
face. As the strength of AF fluctuations is decreased, the 
number of empty energy levels decrease compared to full 
energy levels, the x-band becomes flatter and the Fermi 
surface is bigger (diagram b in Fig. 1). Deep inside the 



heavy Fermi liquid phase, the Kondo effect is dominant, 
and one expects a situation similar to the one impurity 
case, with a totally flat x-band and a large Fermi surface; 
in this phase the x-fermion is no longer a good elemen- 
tary excitation and, x being a flat massless mode, the 
theory is ill defined. Just at the brink of the Fermi liquid 
phase, the x-band is "almost" flat (diagram c in Fig. 1). 
One can thus expect that the velocity vanishes for part 
of the Fermi surface, producing a singularity. 

To summarize, our scenario predicts a new excitation 
close to the QCP, which is fcrmionic in nature and doesn't 
carry a spin. As one integrates out the high energy degree 
of fredom, this new fermion forms a "band": its disper- 
sion has positive and negative energy levels. The band 
structure at the QCP is not universal, it depends on the 
relative strengths of the Kondo eff'ect and AF interactions 
(Fig. 2). If the AF interaction is small compared to the 
Kondo interaction, the QCP lies inside the heavy Fermi 
liquid phase, which is expected to be the case for CeNi2 
Ge2 (diagram a in Fig. 2). In the opposite limit, when 
the AF fluctuations are strong compare to the Kondo ef- 
fect, the Fermi liquid phase can lie outside the AF phase, 
and one can expect singularities at the QCP (diagram 
b in Fig. 2). This is the case for the YbRh2Si2 com- 
pound. The evolution of the x-band structure accounts 
for the variation of the total Fermi surface volume, with 
a smooth cross-over from a rather small Fermi surface in- 
side the AF phase to a big Fermi surface inside the Fermi 
liquid phase. 

We now turn to the speciflc case of YbRh2(Sii_a: Gea;)2 
doped to criticality for x = 0.05. We assume that, at the 
QCP, the x-band is well formed, but the velocity vanishes 
at the hot lines. The effective Lagrangian comprises four 
terms: a conduction electron term Sf, a boson term Sb 
which describes the critical spinous, the x-fcrmion term 
S^, and the interaction Sint between those modes. Both 
the / and x- fermion 's Fermi surface have hot lines re- 
lated by the ordering wave vector Q* (Fig.3). 



S = Sf + S-y^ + Sb + Sint 
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FIG. 3. The /-fermion and x-fermion Fermi surfaces. The 
hot hnes are points related by the ordering wave vector Q*; 
their width is proportional to T^^^ . 

Stnt = 9 j duji duJ2 dPki d^k2 d^q 5q+ki-k2 

where the /-fermion as well as the x-fermion dispersions 
have been linearized around the Fermi surface, the Fermi 
velocity vanishes at the hot lines, leading to a singular 
in k"^ dispersion. We take — a;log(|a;|). It is assumed 
the width of the ^-fermion hot lines varies with tem- 
perature like Aq = {T/aY/^. The spinon dispersion is 
linearized close to the AF QCP , after diagonalizing Hh 



in (1) in a mean field approximation: ujq — y -^^ — 
with A the Lagrange parameter implementing the con- 
straint on the SP(N) representation of the spin, and 
Aq — A(sing2; 4- sing^ -I- sing^). At the QCP, A = A 
and we have linearized LOq around q — (7r/2, 7r/2, 7r/2). 
We first evaluate the dynamical vertex ^(e, O) around 
the hot lines, where e and SI are the incoming fermion 
and boson frequencies respectively and find 



(4) 



The effective theory (3) is thus stable with respect to 
the interaction g. It lies above its upper critical dimen- 
sion. Critical exponents can be evaluated by computing 
self-consistently the self-energies at one loop (Fig. 4) , the 
smallness of the vertex ensuring the validity of this treat- 
ment. Solving the set of three Dyson equations for the 
three fields 6, / and x leads to momentum independent 
self- energies 
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~ i\uj\'^/^ sgn{uj) being irrelevant compared to the 
original damping. 

In order to discuss the thermodynamics of our model, 
we rely on a scaling argument. The free energy is com- 
prised of three terms F = Ff + F^ + Fbos , where Ff is the 




FIG. 4. Feynman diagrams associated with a) the spinon 
self-energy b) the /-fermion self energy c) the x-fermion 
self-energy. Solid, dashed and wavy lines are respectively 
/-fermion, x-fermion and boson propagators. 

free energy of the heavy /-fermions, F^ is the contribu- 
tion from the x fermions and Fi,os from the boson modes. 
Any of these individual contributions is extensive and 
proportional to the phase space volume F ^ Tj^jAg^, 
where Aqi is the momentum width in the i-direction. For 
the x-fermions, only two directions, perpendicular to the 
hot lines, reduce the phase space with Aqi 2 ~ {T /aY^^ 
leading to F^ ^ T{T/aY/^. The boson modes scale in all 
directions, leading to Fbos ~ T{Y.b[T)/Xf ~ T^. The /- 
fermions scale in two directions, perpendicular to the hot 
lines, with Aqi - {T/aY'^ and Aq2 - Y,f{T)/vF, lead- 
ing to Ff,hot ^ T^- Notice that the scaling exponents of 
the hot fermion's and the boson's entropy are the same, 
which is a property of spin-fermion models above their 
upper critical dimension [27]. At low temperatures, F^ 
is the dominant contribution, leading to a specific heat 
coefficient which captures the upturn observed in the ex- 
periments: 



7 = C/T - a^/^r-i/a 



(6) 



We now turn to the behavior of the conductivity at 
criticality. The /-fermions are the sole contributors to 
electric transport because the velocity of the ^-fermion 
vanishes on the hot lines. Assuming that the scattering 
processes from magnetic spinous dominate the transport, 
one distinguishes the contributions from the cold and hot 
regions a{T) = acoid{T) + ahot{T). Following [21], the 
conductivity is the sum of inverse scattering rates asso- 
ciated with the hot and cold regions, multiplied by the 
phase space allowed to these regions 
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where x is a dimensionless parameter which character- 
izes the residual resistivity- for example, the inverse of 
the Residual Resistivity Ratio (RRR), being the ratio of 
the resistivity at 100 K to the resistivity at 10 K. f is a di- 
mensionless parameter characteristic of the temperature 
-one can take for example t = T/Tk, where Tk is the 
one impurity Kondo temperature of the compound, t^^^ 
is the width of the hot lines, is the inverse scattering 
rate of the cold conduction electrons, and t^/^ is the in- 
verse scattering rate of the conduction electrons - up to 
logarithms- from the paramagnetic spinous around the 
hot lines. One sees three regimes in the conductivity: 



t < . ^3/2 ^ ^ ^ ^1/2 . ^1/2 ^ ^ 

p = X + t; cross-over ; Ap ^ . 



(7) 
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FIG. 5. Linear in T resistivity and variation of the tem- 
perature exponent as function oi t = T/Tk- x is the inverse 
Residual Resistivity Ratio. 



With .T ~ 0.1 and ~ 30K for YbRh2Si2, one sees 
that the resistivity is hnear in T, below T = down 
to the lowest temperature (Fig.5), while the behavior 
would appear roughly above T = lOK , corresponding to 
the experimental observation. To conclude the study of 
thermal observables, we focus on the Griineysen param- 
eter, defined as F = a/Cp with a ^ and Cp ^ 
p is the hydrostatic pressure. Calling r = {p—Pc)/Pc the 
departure from the QCP, and noticing, from a small dis- 



placement of the x-fermion hot line that r ~ qfy^ ~ T^/^, 
one finds T ~ ^-2/3 close to the QCP. This exponent is 
in good agreement with the experiments [13] . 

Several theoretical descriptions have been advanced to 
account for the striking experimental results around the 
QCP of heavy fermions. Anisotropic (2D) spin density 
wave scenario [22,23] can explain the linear in T resistiv- 
ity, as well as the logarithmic dependence of the specific 
heat coefficient. A theory with a "local" mode, as well 
as 2D spin fluctuations at criticality, has been advocated 
[23], which also captures the anomalous exponent of the 
spin susceptibility. None of those theories can account for 
the upturn in the specific heat coefficient. Our approach 
has some analogies with the idea of "fractionalization" 
[24] and deconfinement [1], in the sense that a new exci- 
tation appears at the QCP, which is fermionic in nature, 
and which, at high energies, is the gauge field which char- 
acterizes the Kondo bound state. When the new excita- 
tion is fermionic in nature, one might expect the emer- 
gence of a new "phase" where the excitation is stabilized 
[17,25]. Some analogies exist with the "two-fluid model" 
advanced recently to describe CeCoIns [26] : in our model 
also, the impurity spins are partially screened, the differ- 
ence being that at very low energies, the unscreened part 
of the spin, instead of remaining intact, "fractionalizes" 
into a spinon, a spinless fermion, and a "light" electron. 

To summarize, starting from the Kondo-Heisenberg 
model at high energies, we have introduced a new effec- 
tive theory for the QCP of YbRh2(Sio.95 Geo. 05)2- This 
eff'ective theory shows the emergence of a new excitation 
at low energies, fermionic and spinless in nature, charac- 
terizing the Kondo bound state. This new fermion -the 
X-fermion- forms a "band" when high energy degrees of 
freedom are integrated out. The total Fermi surface vol- 



ume evolves from a small Fermi surface to a big Fermi 
surface as the x-fermion band fills up. Assuming a van- 
ishing velocity in the x-fermion band structure close to 
the QCP, one can reproduce many striking experimen- 
tal observations of YbRh2(Sio. 95 Geo. 05)2- One captures 
the upturn in the specific heat coeSicient, with the right 
exponent with temperature, as well as the linear resistiv- 
ity at very low temperatures. The upturn in the specific 
heat coefficient doesn't couple to transport, and couples 
only indirectly -via the spinous and /-fermions- to the 
bulk spin susceptibility. This explains the observation 
that the ratio A/x^ is constant in the heavy Fermi liq- 
uid phase [15]. One also reproduces the variation of the 
Griineysen parameter with the anomalous temperature 
exponent. 
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